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Abstract 

An iterative method is presented to solve the internal and external 
camera calibration parameters, given model target points and their 
images from one or more camera locations. The direct linear transform 
formulation was used to obtain a guess for the iterative method, and 
herein lies one of the strengths of the present method. In all test cases, 
the method converged to the correct solution. In general, an over- 
determined system of nonlinear equations is solved in the least-squares 
sense. The iterative method presented is based on Newton-Raphson 
for solving systems of nonlinear algebraic equations. The Jacobian is 
analytically derived and the pseudo- inverse of the Jacobian is obtained 
by singular value decomposition. 
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1 Introduction 

Certain experimental flow visualization techniques such as the Pressure Sen- 
sitive Paint (PSP), Temperature Sensitive Paint and others involve taking 
images of aircraft models in a wind tunnel. Typically, these images are 
taken with one or more CCD cameras from different positions and orienta- 
tions. These images are then subjected to image registration algorithms and 
mapped onto the model geometry. Details about registration of PSP images 
are discussed by Bell [2]. The mapping of the images to the model geome- 
try, essentially a mapping relationship between model and image coordinate 
systems, is expressed by the equation of photogrammetry (See Marzan and 
Karara [6]). The relationship between the model and image coordinates re- 
quires the determination of the camera calibration parameters. The camera 
calibration paramters are: (a) external or extrinsic parameters which depend 
upon the location and orientation of the camera, and (b) internal or intrinsic 
parameters which are the camera focal length, other lens parameters, dis- 
tortion parameters, etc. The camera parameters which are relevant to this 
report are discussed in more detail in Section 2. Typically, certain control 
points on the model and their projected images are known a priori. Given the 
coordinates of the control points, the equations of projective photogramme- 
try, which are nonlinear equations, are solved to yield the camera orientation. 
The most common method to solve the equations with an iterative solution 
technique is termed bundle adjustment ([5, 4]). Another technique to ob- 
tain the mapping between the image and model space is the Direct Linear 
Transform (DLT) method [1]. 

In this report, a method to solve for the camera calibration parameters 
is described. The solution method uses the DLT to obtain a guess for the 
camera parameters, and then iteratively solves the nonlinear equations us- 
ing a Newton method. The Jacobian for the Newton iterative technique is 
obtained analytically, and the inverse (or rather pseudo-inverse) Jacobian is 
obtained by singular value decomposition. 


4 


2 Solution Method 

2.1 Equations of Photogrammetry 

Let the model and image coordinate systems be denoted by (X,Y,Z) and 
(x, y, z), respectively. The origin of the image coordinate system is located 
at ( Xq , Yci Z c ) model coordinate space. The orientation of the image 

coordinate system with respect to the model system is given by a 3 x 3 
rotation matrix R. The nine terms i =1,2,3 j = 1,2,3 of R are 
not independent and can be reduced to three parameters by taking the Euler 
angle approach, i.e., R is a combination of three separate rotations about the 
three principal axes. Let the rotation about the x, y and z axes be denoted 
by u>. and k, respectively. Then 

R= R z (K)Ry(ct>)Rx{u) (1) 

where each individual rotation matrix is given by 

'10 O' 

R x (oj) = 0 cosw — sin a; , Ry{(p) 

0 sin u cos w 

cos n — sin k 0 
R z (k) — sin k cos rc 0 
0 0 1 

One must be careful and note that the rotations do not commute. After 
matrix multiplication, R is given by 

" cos (p cos k — cos u) sin k + sin u) sin 4> cos k sin u> sin k + cos w sin (p cos k, 
R = cos (p sin k cos u cos k + sin u sin <p sin k — sin uj cos k + cos to sin (p sin k 

— sirup sin u cos (p cos w cos (p 

Note that the rows (and columns) of the above matrix form orthogonal bases 
vectors, a property which will be exploited later. 

In actual practice, several target locations are identified on the model and 
their coordinates are measured. The model is imaged by a CCD camera. For 
an ideal camera, the projection of a model point (X, Y, Z ) is ( x , y) in image 
coordinates. However, due to image distortion, the image of the model point 
is measured at (x\y'). 


cos <p 0 sin cp 
0 1 0 
— sin <p 0 cos (p 
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The relationship between model and image coordinates is expressed by 
the projective equation of photogrammetry[6] 

x — x p = x' — x p + Ax 

_ , T n(A~ — Xg) + £ 12 ~ Yp) + r 13(Z ~ ^c) 

^ 3 i (A r — Xq) + r 32 (Y ~ Yc) + r 33 (Z — Zp ) 
y — y p = y' — y P + Ay 

_ , 7*21 (AT — Xg) + r 22 (Y — lc) + r 2 3 (X — Zc) 

r 31 (X - X c ) + r 32 (Y - *c) + r 33 (Z - Z c ) - 

In equation (2) Ax, Ay are the terms which model the effects of symmet- 
rical and asymmetrical lens distortion; x p ,y p are the image coordinates of 
the principal point; and rp- is the i th row and j th column element of R ; 
and / is the focal length of the camera. A commonly used set of functions 
for Ax and Ay is 

Ax = x{kir 2 + k 2 r 4 + /c 3 r 6 ) + p x (r 2 + 2x 2 ) + 2 p 2 xy. 

Ay = y(Aqr 2 + k 2 r 4 + k 3 r 6 ) + 2pixy + p 2 {r 2 + 2 y 2 ) + a x x + a 2 y, (3) 

where x = x' - x p , y = y' - y p , r 2 = x 2 + y 2 ; and ki, Pi and a x rep- 
resent radial, lens decentering, and affinity distortion terms, respectively. 
The six-tuple (X c , Y c , Z c , w, <f>, k) which gives the position and orientation 
of the camera contains the external calibration parameters ; while the ten- 
tuple {ki,k 2 ,k 3 ,p u p 2 , ai,a 2 ,x p ,y p ,f) representing the radial, lens decenter- 
ing, affinity distortion terms, the principal point, and the focal length con- 
tains the internal calibration parameters. 

Let the number of model target points be m. A given camera produces 
images of the model from n different locations. For each camera location, the 
image coordinates of a subset of the m, model points may be measured. Let 
the number of measured image points for the i th camera location be mi < m. 
The total number of unknowns are 6n + 10 where 6n are the six external 
calibration parameters for n camera locations and 10 internal calibration 
parameters. The total number of equations available is 2 £^=1 m i < 2nm. 
For all the test cases discussed later in this report, we encounter an over- 
determined system, i.e., 2£" =1 nq > 6n + 10. 

2.2 Solution Method 

The over-determined system of equations governing the internal and external 
calibration parameters is solved in a least-squares sense as follows. We first 
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recast the equations in a convenient form. Define a vector of unknowns as 

£ = Z(Xc,u * 1 , Ar c>2 , $ 2 , • • • Xcjn, ^) T (4) 

where X c ,i = (X c ,uY Ciii Z c ,i) and <1\ = {u) U <f> u kJ are the camera external 
calibration parameters for the i th camera location; and 

= (k u k 2 ,k 3 ,p u p 2 ,aua 2 ,x p ,y p ,f) is the vector of the camera internal 
calibration parameters. Then rewrite the equations as 

T 7 (C\ - f F ‘ (0 = + " *!>)(** " *c,i) • h,i + f{Xj - X c ,i) ■ r h i = 0 

^ ~ [ G,(0 = (</' + - a: p ) ( Aj - X Cli ) • f 3 , t + f(Xj - X C J ■ f M = 0 ( jjj 

In equation (5), l = + j, (mo = 0); j — 1, 2, ■ ■ -m*; and i — 1, 2, ■ • -n, 

with the index j varying faster than index i. Furthermore, Aj is the position 
vector of the j th target point, i.e., Xj = ( Xj,Yj , Zj ); and X c ,i is the position 
vector of the i th camera, i.e., X c ,i = (X c ,i,Y Ct i, Z c , i)- The unit vectors 
f fc>i correspond to the k th column vector of the rotation matrix R and are 
functions of the i th camera Euler angles (w*,^, k,,). 

The above system of nonlinear equations is solved in an iterative fashion 
as described below. Assume the solution vector £ is known at iteration ’ k\ 
The left-hand side of equation (5) may be expressed as a truncated Taylor 
series at iteration k + 1. 

^* +1 (0 = ^*(0+(^) A£*, (6) 

where A£ fc = £* +1 - £ fc . In the above equation, T is a column vector with 
2 Yl™ m i elements, the Jacobian dT/d^ is a matrix with 2^"mj rows and 
(6n + 10) columns, and A£ is a column vector with (6n + 10) elements. 
Following the Newton-Raphson iterative technique, we set T k+1 = 0 and 
solve for £ k+l as 

e*'=e- /*(«). m 

where the term [] -1 is the inverse (pseudo-inverse) of the Jacobian matrix 
if 2E ?mi = 6 n + 10, (2 > On + 10). The elements of the Jacobian 

matrix are analytically determined (see Appendix A for details). At each 
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iteration, we solve a system of linear equations. For over-determined systems 
(2 m l > 6n+ 10) the linear system is solved in the least-squares sense. The 
above iterative procedure is terminated upon convergence. It is well known 
that the standard Newton-Raphson method converges if the initial guess is 
sufficiently close to the solution. The initial guess for the above procedure is 
discussed in detail in section 2.3. 


2.2.1 Pseudo-inverse and Singular Value Decomposition 

For an over-determined system (2 Yh m i >6 n + 10), the Jacobian matrix 
is not square. A preferred method for solving linear least-squares problems 
is singular value decomposition (SVD)[9]. In general, any M x N matrix A 
may be written as the product of an M x N column-orthogonal matrix U, 
an N x N diagonal matrix W, and the transpose of an N x N orthogonal 
matrix V. The system of equations A.x = b is solved in the least-squares 
sense by finding x which minimizes x 2 = ||A ■ x — 6||. The solution x is given 
by 


x = V- [diagil/wj)} - U T -b , (8) 


where Wj is the jth diagonal element of the diagonal matrix If and are called 
the singular values of A. The product of the three matrices, 

V ■ [ diag(l/wj )] • U T , in the equation above is the pseudo-inverse of A. In 
our case, the pseudo-inverse of the Jacobian is calculated using SVD, and 
the solution vector at £* +1 is calculated in the least-squares sense as outlined 
above. An added benefit of using singular value decomposition is that the 
above procedure works even when the Jacobian is singular or ill-conditioned. 
We compute the condition number (defined as the ratio of the largest to the 
smallest singular value) of the Jacobian matrix. If the inverse of the condition 
number is smaller than 10~ p , all singular values which are smaller than 10 _p 
times the maximum singular value are set to zero. This prevents round-off 
error from contaminating the solution. In our case, we choose p = 12 because 
our computations are done using double precision arithmetic. 


2.2.2 Convergence Criteria 

The convergence criteria to terminate the iterative solution is 

|^rfc+l _ jr*| < 6j 


( 9 ) 
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where e is a “small” number determined by the user. Another alternative is 
to examine the change in the norm of the solution vector For both 

criteria, convergence is achieved in 10 — 20 iterations for e = 10“ 8 . In no 
test case did we need to iterate more than 20 times. Convergence was not 
quadratic, as one would expect for a standard Newton- Raphson technique. 
Perhaps this is not surprising because we are solving the equations in a 
least-squares sense, and in some cases the Jacobian matrix may have large 
condition numbers. 


2.2.3 Gimbal Lock 


We now focus attention on the case when $ = 7t/2. The rotation matrix R 
for <j) = 7t/2 becomes 


R = 


0 

0 

-1 


sin (a; — k) 
cos (a; — k) 
0 


cos(w — k) 
— sin (a; — n) 
0 


( 10 ) 


This shows that the solution only depends upon the difference of ui and k, 
i.e., the number of independent unknowns reduces by one. This is called 
“Gimbal Lock”. Further discussion of this issue is- given in Section 2.4. 


2.3 Initial Guess 

For a good optical system that does not suffer from optical distortion, the 
error terms in equation (2) may be neglected. Some algebraic manipula- 
tion leads to the following form, called the direct linear transform (DLT) 
formulation [1], 

L x X + L 2 Y + L 3 Z + L 4 L5X + LeT + L7Z + L 8 qjx 

X LqX + LiqY + L\\Z ' ^ LqX + L\qY + L\\Z 

where Li---L n are called the DLT coefficients. Given the model target 
points and their images, one solves for the DLT coefficients using the solution 
procedure given in Reference [8]. The DLT coefficients for n camera locations 
are obtained. 

2.3.1 Estimate of Camera Position 

In equation (2) the numerator and denominator are zero w r hen X = X c , Y = 
Y c , Z = Z c . This implies that the camera position vector X c may be deter- 
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mined by solving the following linear system of equations: 

LiXq + L 2 Yq + L%Zc + = 0, 

L$ Xc + L§Yc + LjZq + L$ = 0, 

L$Xc + LiqYq + L\\Zq + 1 = 0. 


(12) 


2.3.2 Estimate of Camera Focal Length 


Recall that one of the internal calibration parameters to be determined is 
the camera focal length /. For convenience we define some new variables as 
follows: 


, _ L \ _ ,Di , _ L 2 _ Tn 

sl — T J js2 — r I 3 

T 31 L g r 31 


T)i = - — = , 

Mo r 32 


^2 ,H2 

V 2 = ~j— = -/ — > 

^10 r 32 


n2 


r= — = -f 

>1 r ^ r ’ r * rnn 

Ln r 33 ^33 


c — ^3 _ r r 13 

S3 r J 1 

L>9 T 31 

£3 r r 13 

% = T J i 

±>10 r 32 

> _ ^3 _ r r 13 

S3 = 7 — — — / ■ 

^>11 r 33 


(13) 


Now we use the fact that the row and column vectors of R are orthogonal 
unit vectors. It can then be trivially shown that 


J_ - i_ _L _L 
7 5- I%? + < V 


(14) 


where = £ 2 + and so on. The focal length is then calculated 

by taking the positive square root of / 2 . In the estimation of the focal 
length, an alternative is to use L 4 , L 5 , L 6 and r 2 1 , r 2 2 , r 33 instead of Li,L 2 , L 3 , 
r u , r 12 , r 13 , respectively, in equation (13). Finally, because we have n camera 
locations, we calculate n sets of DLT coefficients and n focal lengths. We 
choose the starting point of the focal length of the camera as the average of 
the n estimated focal lengths. 


2.3.3 Estimate of Camera Orientation 

We begin by using the last row of R, the calculated DLT coefficients, and 
the estimated focal length to estimate the orientation Euler angles (u, q k ) 
for each camera location. We first estimate <j> by the following relation 

sin 2 (j) = 4 


(15) 
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We take the positive square root and take the principal value for 0. The 
other two angles are estimated by using the following equations: 


sin 2 cj = 


P 

rtf cos 2 4 > 1 


COSK - 


Li sin (p 

JLg COS ( f )' 


(16) 


For the rare case of 0 = 7t/2> we simply set (j) = 0.957t/2. This is justified 
because we are simply interested in an initial guess for <fi. This eliminates 
the possibility that the denominator in equation (16) is zero. Note that 
as calculated above (j) £ [0,7r/2], u £ [0, vr/2] and k £ [0,7r]. Using these 
calculated angles, one may recalculate the focal length using the following 
relations: 


/ = 


L i7 

T-ii 


/ = 


L 57 

7-21 ’ 


, _ -^27 f-PPL f = ~^47 

J T \2 ’ r 13 T\\Xc + t u Yc + T\$Zc 

,_Vy t-hl f = 

T 22 C23 ’ + r n Yc + r 23^c' 


(17) 


where 7 = r^Xc + + ^33 Zc- The estimate of the angles are refined 

by using —9, n ±6,9 - it, (9 = a 0, /c), to calculate R which is then used 
in equation (17). That combination of angles which gives a positive focal 
length for all the above eight calculations in equation (17) is taken as the 
initial guess of the orientation. 


2.4 Test Cases and Uniqueness of Solution 

For nonzero Ax and Ay in equation (2), there are a multiplicity of solutions. 
It is unclear how many solutions are actually present. Contrived tests were 
carried out with known camera calibration parameters. When the above 
procedure to obtain an initial guess of the solution is adopted, the solution 
obtained was always the correct one. For some contrived tests, if the so- 
lution vector was initialized to arbitrary quantities, the converged solution 
was not the correct one. In some cases, the solution converged to a very 
small focal length with large distortion parameters. It is recommended that 
the image points be recomputed once the camera calibration parameters are 
determined. Then a calculation of the root-mean-square error between the 
computed and the given image point coordinates is performed. For model 
and target points provided by Bell [3] for the B-70 and Boeing airplane test 
data, the root-mean-sqnare error w r as always less than one pixel. 
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If Ax = Ay = 0 and if the focal length is provided, then the above 
iterative procedure converges to a unique solution. Of course, the solution 
for the orientation angles is unique up to an additive constant. It is trivial 
to ensure that the angles always lie between — 2n and 2n. 

For cases of Ginibal lock, tests indicated that the solution converges to 
arbitrary values of u> and /c, while the correct difference between these two 
angles is always attained. Thus, nothing special is done for the case when 0 = 
7 t/ 2. The reader is simply cautioned that the difference uj—k is important and 
not the individual values. Because we are eventually interested in mapping 
between the image and model spaces, the mapping obtained even in the case 
when we have gimbal lock is correct. 

Finally, we note that the above solution procedure was coded into a C- 
program library which may be obtained by sending electronic mail to the 
author. The library was made flexible enough so that, in the event the 
code converges to an undesirable solution (which is discerned quickly by an 
examination of the root-mean-square error between measured and calculated 
image point locations), the user may specify an initial guess for the solution. 
We now present two specific sample calculations. The data used in these 
sample calculations is presented in tabular form, and the interested reader 
may obtain the numerical values by sending electronic mail to the author. 

2.4.1 Example 1: Joukowski Wing Dataset 

A dataset was artificially generated to mimic the situation encountered in 
pressure sensitive paint experiments. Basically, we created a wing using sev- 
eral two-dimensional cross sections, each of which was a Joukowski airfoil [7]. 
Then we sprinkled 46 target control points in a random fashion on this wing. 
A camera was placed at four different stations. The resulting four images, 
each of size 512 x 512 pixels, are shown in Figure 1. The target points, and 
their locations in the images are given in Table (1), and Tables (2)-(5), re- 
spectively. Note that the image coordinates, in pixels, are usually scaled by 
a certain factor in order to operate in the same units as the model. In this 
example, the scale factor was 1/512. The above method was used to deter- 
mine the camera calibration parameters and the results are given in Table 
(6). In this contrived example, in 7 iterations \T k+l — T k \ = 5.7 x 10 8 and 
in 8 iterations |^ fc+1 - F k \ = 1.3 x 10~ 15 . Furthermore, in this example in 
8 iterations we also achieved |.F* +1 | < 10~ 15 . Another point to note in this 
example, is that, for station numbered T, the orientation angle (j) = 90° and 
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therefore we have a gimbal lock. The images were generated using u = 0° 
and k = 50°. The solution obtained shows u = 1.381° and k = 51.381° so 
that the difference between these angles is still 50° as it should be for the 
gimbal lock case. The internal parameters are very small as they should be 
in this example. 

2.4.2 Example 2: Boeing Airplane Dataset 

Data in the form of target point locations and their images was provided by 
Bell [3] for a Boeing airplane. The target points are given in Table (7). Of 
course, the units are normalized so that they do not correspond to the real 
dimensions of the actual model in the experiment. Three images, each of size 
512 x 1024 pixels, were provided. The coordinates of the target points for each 
of the three images are given in Tables (8)-(10). The image coordinates were 
scaled by a factor 9.4488 x 10 -4 . The above method was used to determine 
the camera calibration parameters, and the results are given in Table (11). 
In this example, convergence is achieved in 10 iterations w'hen we obtain 
|JF* +1 — T k | = 3.66 x 10 -9 . The L norm of T settled down to 0.051818 after 
6 iterations. Also shown in Tables (8)-(10) are the calculated coordinates of 
the target points in the images using equation (2). The root-mean-square 
error between the calculated and the measured locations of the image points 
is less than 0.2 pixels. 


3 Conclusion 

In this report, we have presented an iterative method to solve the internal 
and external camera calibration parameters given model target points and 
their images from one or more camera locations. The direct linear transform 
formulation was used to obtain a guess for the iterative method, and herein 
lies one of the strengths of the present method. In all test cases, the method 
gave the correct solution without user intervention in selecting a good initial 
guess. Sufficient details are given In the report, including detailed analytical 
expressions for the Jacobian, to enable one to write a computer program. 
Alternatively, a C-library which uses the present method to determine the 
camera calibration parameters may be obtained by sending electronic mail 
to the author. 



Figure 1 : Images corresponding to the four different camera locations. The 
black dots on the images correspond to the target points. The color corre- 
sponds to the coefficient of pressure for the wing with red (blue) correspond- 
ing to high (low) values. 
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4 Appendix A: Derivation of the Jacobian 


The Jacobian matrix in equation (6) is given by 


dT 

d£ 


dF 

dXni 

dd 


dX C ,\ 

0 

0 

0 

0 


dF 

a<s>! 

dG 

d$i 

0 

0 


0 

0 


0 0 

0 0 

dF dF 


dXc,2 d$2 
dG dG 


dXca d*f>2 

0 0 

0 0 


0 

0 

0 

0 


dF 

dx Cn 

ad 

dX C>n 


0 

0 

0 

0 


d£ "I 

dv 

dG 

d<b 

dF 

dV 

dG 

d ¥ J 


dF dF 

a * 

dG dG 

d$n a* 


(18) 


Each row of the Jacobian matrix consists of 6n + 10 entries, of which 6 n — 6 
entries are zero. The derivatives of F (and G ) with respect to Xcj, 4>j and 
T are row vectors with 3, 3, and 10 elements, respectively. These derivatives 
are given below. 


dF 

dX c 

dG 

dX c 


- -(x - x p )f 3 - ffi, 


= ~{y - y P )h - ffi, 


(19) 


OF 

5<f> 

dG 

5<f> 


(x-x,)(X-Xc)-^ + f(X-X c )^, 
(y-y,)(x-x c )-^ + i(x-Xc)% 


( 20 ) 


where 




dr i 

0 

sin cd sin k + cos to sin 0 cos k 
cos c o sin k — sin u sin 0 cos k 

— sin 0 cos k 
sin u) cos 0 cos k 
COS U) cos 0 COS K 

— cos 0 sin k 

— cos u) cos k — sin u sin 0 sin k 
sin u cos k — cos oo sin 0 sin k 
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dr 2 


0 — sin (f) sin k cos 4> cos k 

— sin u cos k + cos u sin <fi sin k sin u cos (p sin k — cos u sin k + sin u sin (p cos k 

— cos u cos k — sin w sin <f> sin k cos w cos <p sin k sin u sin k + cos ui sin (p cos k 


df 3 

<9$ 


0 — cos <p 0 

cos lo cos (p — sin u sin <p 0 
cos lo cos <p sin k — cos ui sin (p 0 


The derivatives with respect to the internal calibration parameters are given 
below: 


— = (X - X c ) • r 3 xr , — = (X - X c ) ■ h xr , — = (X - X c ) ■ r 3 xr b , 


dF 

dk { 

dF 

dp i 

dF_ 
da i 
dF 
dx p 
dF_ 
dy P 
dF 
df 

and 

dG 

dki 

dG 

dp2 

3G 

da i 

dG 


0x p 

dG_ 

dy P 

dG 

df 


dF 


dF 


dh 


dk 3 


dF 


{X - X c ) ■ h ( r 2 + 2a: 2 ), ~ = {X - X c ) ■ f 3 2 xy, 

dF 

0, ^- = o, 

0(12 


= -{X - X c ) • f 3 1 


= ~{X — X c ) • r 3 


— = {X-Xc)-f u 


dAx\ 

dx p J ’ 


3Aa: 
dy P ’ 


(21) 


3C BC 

= (X - x c ) ■ r 3 yr 2 , — = (X - X c ) ■ r 3 W 4 , — = (<T - ^c) • 

BC 

= (X - X c ) ■ r 3 (r 2 + 2?J 2 ), — — = (X - X c ) • r 3 2xy, 


dpi 


dG 


(X - X c ) hx, — = (X - X c ) ■ r 3 y, 


_ -(X-Xc)-f 3 


dAy 
dx p ’ 




^ = (X-X c ).f 2 


( 22 ) 
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(hi + k 2 r 2 -I- k 3 r 4 )r 2 - 2x 2 (ki + 2 k 2 r 2 + 3k 3 r 4 ) - 6pirr - 2 p 2 y, 
2xy(k\ + 2k 2 r 2 + 3k 3 r 4 ) — 2p\y — 2p 2 x, (23) 

2xy{k\ + 2 k 2 r 2 + 3 k 3 r 4 ) - 2 p x y - 2 p 2 x - a x , 

(ki + k 2 r 2 + k 3 r 4 )r 2 - 2y 2 (k 1 + 2 k 2 r 2 + 3 k 3 r 4 ) - 2 p x x - Qp 2 y - a 2 . 
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Table 1: Target locations for the Joukowski wing dataset (Example 1). 


No. 

Target Locations 

No. 

Target Locations 


X 

Y 

Z 


X 

X 

Z 

1 

3.205 

-0.097 

-12.750 

24 

-2.144 

0.062 

-6.750 

2 

2.289 

-0.107 

-3.375 

25 

-2.054 

0.124 

-9.000 

3 

2.099 

-0.135 

-1.875 

26 

-1.899 

0.224 

-11.625 

4 

2.134 

-0.225 

-3.375 

27 

-1.614 

0.332 

-12.375 

5 

2.925 

-0.651 

-14.250 

28 

-1.121 

0.211 

-1.875 

6 

2.397 

-0.587 

-10.125 

29 

-0.897 

0.352 

-5.250 

7 

1.505 

-0.190 

-0.750 

30 

-0.587 

0.586 

-10.125 

8 

1.829 

-0.589 

-8.250 

31 

-0.197 

0.824 

-13.875 

9 

1.225 

-0.291 

-2.625 

32 

0.192 

0.951 

-14.625 

10 

1.095 

-0.386 

-4.500 

33 

0.339 

0.711 

-7.875 

11 

0.648 

-0.162 

-1.125 

34 

0.674 

0.789 

-8.625 

12 

0.910 

-0.781 

-13.125 

35 

1.124 

0.966 

-11.625 

13 

-0.304 

-0.334 

-8.250 

36 

1.280 

0.840 

-9.000 

14 

-1.007 

-0.315 

-12.750 

37 

1.030 

0.416 

-1.125 

15 

-1.356 

-0.170 

-10.500 

38 

2.097 

0.990 

-13.125 

16 

-1.558 

-0.022 

-4.875 

39 

2.090 

0.779 

-9.750 

17 

-2.038 

-0.062 

-11.625 

40 

2.391 

0.754 

-10.875 

18 

-1.827 

0.019 

-1.125 

41 

2.839 

0.766 

-13.875 

19 

-2.035 

0.011 

-2.625 

42 

2.131 

0.352 

-4.125 

20 

-2.222 

0.005 

-4.125 

43 

1.939 

0.185 

-0.750 

21 

-2.098 

0.001 

-1.500 

44 

3.168 

0.348 

-13.125 

22 

-2.354 

0.001 

-5.250 

45 

2.795 

0.162 

-8.625 

23 

-2.454 

0.009 

-7.500 

46 

2.999 

0.048 

-10.500 




19 


Table 2: Coordinates of target points in the image collected by the camera 
at station 1 for the Joukowski wing dataset (Example 1). The image is of 
size 512 x 512, and the image coordinates are in pixels. 


No. 

Image Coordinates 

No. 

Image Coordinates 


X 

y 

X 

y 

1 


146.457 

24 

160.456 

224.876 

2 



25 

176.300 

204.277 

3 

95.147 


26 

194.569 

179.656 

4 



27 

198.479 

171.013 

5 

212.553 

136.367 

28 

116.458 

267.469 

6 

174.584 


29 

140.216 

234.806 

7 


284.589 

30 

175.208 

186.178 

8 


207.777 

31 

201.894 

146.481 

9 


265.472 

32 

205.808 

135.855 

■a 


246.778 

33 

151.325 

204.287 

ii 


279.287 

34 

155.103 

195.114 

12 



35 

177.012 

160.952 

13 

167.772 


36 

154.601 

188.768 

14 



37 

92.015 

274.746 

15 

187.955 


38 

186.246 

139.137 

16 


241.727 

39 

157.572 

178.074 

17 

197.497 

182.269 

40 

166.679 

164.445 

18 


275.284 

41 

192.987 

127.665 

19 


261.814 

42 

110.244 

243.567 

wm 

141.855 

248.571 

43 

82.843 

281.860 

21 



44 

189.287 

137.977 

22 


238.776 

45 

149.809 

193.604 

23 

167.733 

219.143 

46 

168.010 

172.339 
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Table 3: Coordinates of target points in the image collected by the camera 
at station 2 for the Joukowski wing dataset (Example 1). The image is of 
size 512 x 512, and the image coordinates are in pixels. 


No. 

Image Coordinates 

No. 

Image Coordinates 


X 

y 


X 

y 

1 

273.318 

266.992 

24 

196.820 

341.550 i 

2 

253.257 

393.849 

25 

201.420 

310.587 

3 

249.113 

416.495 

26 

207.185 

276.344 

4 

251.268 

394.358 

27 

211.936 

267.293 

5 

272.096 

246.787 

28 

204.016 

412.877 

6 

261.764 

299.804 

29 

212.029 

362.831 

7 

239.547 

434.061 

30 

222.415 

297.045 

8 

252.447 

324.952 

31 

231.513 

251.418 

9 

237.834 

405.564 

32 

236.943 

243.310 

10 

238.230 

377.750 

33 

231.682 

327.042 

11 

227.792 

427.650 

34 

236.750 

317.417 

12 

245.354 

258.497 

35 

245.249 

280.090 

13 

223.538 

322.788 

36 

244.725 

313.051 

14 

219.798 

262.095 

37 

233.126 

423.883 

15 

212.420 

291.232 

38 

258.276 

263.039 

16 

201.994 

369.081 

39 

255.563 

304.222 

17 

205.038 

275.874 

40 

260.278 

290.383 

18 

192.432 

425.380 

41 

268.050 

254.517 

19 

191.853 

402.181 

42 

251.325 

380.749 

20 

191.578 

379.572 

43 

245.391 

431.519 

21 

189.092 

419.533 

44 

272.227 

263.058 

22 

191.506 

363.001 

45 

264.379 

319.792 

23 

193.621 

330.905 

46 

268.712 

295.485 
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Table 4: Coordinates of target points in the image collected by the camera 
at station 3 for the Joukowski wing dataset (Example 1). The image is of 
size 512 x 512, and the image coordinates are in pixels. 


No. 

Image Coordinates 

No. 

Image Coordinates 


X 

y 


x 

y 

1 

337.370 

315.663 

24 

275.158 

235.344 

2 

338.494 

193.709 

25 

274.256 

263.767 

3 

338.016 

172.490 

26 

273.800 

296.027 

4 

336.123 

193.743 

27 

276.530 

305.758 

5 

331.017 

329.827 

28 

293.282 

170.472 

6 

329.749 

280.974 

29 

292.566 

216.949 

7 

331.159 

156.270 

30 

291.402 

280.842 

8 

324.820 

257.256 

31 

292.451 

328.188 

9 

324.480 

182.952 

32 

296.587 

338.494 

10 

320.130 

208.458 

33 

306.162 

253.432 

11 

318.596 

161.344 

34 

309.823 

263.852 

12 

307.405 

313.285 

35 

312.269 

303.571 

13 

297.457 

255.452 

36 

317.657 

269.655 

14 

284.349 

308.228 

37 

324.532 

159.827 

15 

281.972 

281.989 

38 

323.479 

323.971 

16 

284.542 

211.460 

39 

327.750 

280.318 

17 

272.525 

294.754 

40 

330.297 

295.116 

18 

284.647 

160.137 

41 

332.208 

333.185 

19 

280.418 

180.911 

42 

335.878 

203.741 

20 

276.576 

201.061 

43 

338.005 

155.220 

21 

280.674 

165.361 

44 

337.050 

322.356 

22 

273.887 

215.788 

45 

338.187 

264.851 

23 

270.691 

244.470 

46 

338.046 

288.547 
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Table 5: Coordinates of target points in the image collected by the camera 
at station 4 for the Joukowski wing dataset (Example 1). The image is of 
size 512 x 512, and the image coordinates are in pixels. 


No. 

Image Coordinates 

No. 

Image Coordinates 


X 

y 


X 

y 

1 

337.370 

315.663 

24 

275.158 

235.344 

2 

338.494 

193.709 

25 

274.256 

263.767 

3 

338.016 

172.490 

26 

273.800 

296.027 

4 

336.123 

193.743 

27 

276.530 

305.758 

5 

331.017 

329.827 

28 

293.282 

170.472 

6 

329.749 

280.974 

29 

292.566 

216.949 

7 

331.159 

156.270 

30 

291.402 

280.842 

8 

324.820 

257.256 

31 

292.451 

328.188 

1 9 

324.480 

182.952 

32 

296.587 

338.494 

10 

320.130 

208.458 

33 

306.162 

253.432 

11 

318.596 

161.344 

34 

309.823 

263.852 

12 

307.405 

313.285 

35 

312.269 

303.571 

13 

297.457 

255.452 

36 

317.657 

269.655 

14 

284.349 

308.228 

37 

324.532 

159.827 

15 

281.972 

281.989 

38 

323.479 

323.971 

16 

284.542 

211.460 

39 

327.750 

280.318 

17 

272.525 

294.754 

40 

330.297 

295.116 

18 

284.647 

160.137 

41 

332.208 

333.185 

19 

280.418 

180.911 

42 

335.878 

203.741 

20 

276.576 

201.061 

43 

338.005 

155.220 

21 

280.674 

165.361 

44 

337.050 

322.356 

22 

273.887 

215.788 

45 

338.187 

264.851 

23 

270.691 

244.470 

46 

338.046 

288.547 




Table 6: Camera parameters for the Joukowski wing dataset (Example 1). 
Note that the RMS error and the principal locations ( x p ,y p ) are in pixels, 
while the orientation angles are in degrees. 


Parameter 

Station 1 

Station 2 

Station 3 

Station 4 

Xc 

20.000 

0.000 

0.000 

-10.000 

Y c 

-8.000 

-20.000 

20.000 

10.000 

Z c 

-10.000 

-10.000 

-5.000 

0.000 

UJ 

1.381 

-80.000 

80.000 

60.000 

0 

90.000 

5.000 

-10.000 

30.000 

K 

51.381 

5.000 

5.000 

45.000 

ki 

-4.053E-7 

-4.053E-7 

-4.053E-7 

-4.053E-7 

k-2 

-2.179E-6 

-2.179E-6 

-2.179E-6 

-2.179E-6 


4.107E-5 

4.107E-5 

4.107E-5 

4.107E-5 

Pi 

1.052E-9 

1.052E-9 

1.052E-9 

1.052E-9 

P2 

8.859E-8 

8.859E-8 

8.859E-8 

8.859E-8 

a x 

-2.334E-9 

-2.334E-9 

-2.334E-9 

-2.334E-9 

a 2 

-6.064E-9 

-6.064E-9 

-6.064E-9 

-6.064E-9 

Xp 

1.000 

1.000 

1.000 

1.000 

Vp 

1.000 

1.000 

1.000 

1.000 

/ 

0.500 

0.500 

0.500 

0.500 

x RMS Error 

5.040E-6 

3.470E-06 

4.070E-06 

3.535E-06 

y RMS Error 

3.866E-6 

2.788E-06 

2.388E-06 

6.169E-06 
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Table 7: Target locations for a Boeing airplane dataset (example 2). 


No. 

Target locations 


X 

Y 

Z 

1 

1607.182 

115.070 

227.882 

2 

2032.488 

97.430 

231.272 

3 

2034.920 

262.636 

221.354 

4 

2290.974 

93.098 

231.192 

5 

2289.736 

336.418 

219.346 

6 

2456.156 

235.920 

226.210 

7 

2573.540 

92.632 

224.990 

8 

2577.158 

422.288 

220.214 

9 

2820.132 

208.338 

215.554 

10 

2783.998 

373.076 

219.732 

11 

2717.762 

551.874 

221.234 

12 

2874.836 

741.880 

221.912 

13 

1903.644 

245.944 

213.790 

14 

2088.480 

322.896 

197.480 

15 

2178.304 

350.492 

198.340 

16 

2426.106 

435.612 

198.062 

17 

2645.676 

604.052 

202.698 

18 

2869.044 

568.930 

213.628 
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Table 8: Measured and calculated image locations, in pixels, for camera at 
station 1 for the Boeing dataset (Example 2). 


No. 

Measured 

Calculated 


X 

y 

%calc 

Ucalc 

1 

462.362 

103.700 

462.369 

103.701 

2 

472.600 

367.667 

472.791 

367.531 

3 

373.792 

370.304 

373.621 

370.230 

4 

473.892 

521.270 

473.837 

521.459 

5 

331.696 

519.208 

331.713 

519.089 

6 

389.500 

614.500 

389.474 

614.347 

7 

471.167 

680.571 

470.999 

680.585 

8 

285.167 

677.500 

285.035 

677.350 

9 

403.333 

807.333 

403.382 

807.255 

10 

313.056 

786.139 

313.192 

786.158 

11 

216.000 

748.083 

216.262 

748.666 

12 

120.611 

822.389 

120.641 

822.225 

13 

383.222 

291.444 

383.156 

291.502 

14 

339.267 

401.867 

323.947 

454.114 

15 

323.783 

454.317 

339.303 

402.114 

16 

277.933 

592.467 

277.922 

592.604 

17 

189.500 

707.500 

189.416 

707.300 

26 

210.500 

822.500 

210.319 

822.526 
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Table 9: Measured and calculated image locations, in pixels, for camera at 
station 2 for the Boeing dataset (Example 2). 


No. 

Measured 

Calculated 


X 

y 

calc 

Ucalc 

1 

455.778 

154.861 

455.800 

154.794 

2 

471.736 

358.264 

471.947 

358.090 

3 

374.000 

363.583 

373.839 

363.730 

4 

478.204 

501.444 

478.263 

501.570 

5 

327.167 

503.429 

327.124 

503.329 

6 

387.903 

601.097 

387.768 

601.140 

7 

481.600 

677.000 

481.429 

676.947 

8 

267.042 

675.137 

266.967 

675.094 

9 

403.000 

840.905 

403.168 

840.768 

10 

292.892 

811.175 

292.935 

811.069 

11 

179.104 

761.340 

179.548 

761.946 

12 

51.852 

858.571 

51.845 

858.545 

13 

384.111 

299.444 

384.075 

299.653 

14 

338.889 

398.556 

338.954 

398.631 

15 

321.444 

446.444 

321.520 

446.285 

16 

265.056 

586.333 

265.005 

586.267 

17 

153.389 

716.611 

153.220 

716.359 

18 

162.944 

861.944 

162.700 

861.915 
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Table 10: Measured and calculated image locations, in pixels, for camera at 
station 3 for the Boeing dataset (Example 2). 


No. 

Measured 

Calculated 


X 

y 

%calc 

Vcalc 

1 

456.611 

144.778 

456.536 

144.796 

2 

471.905 

359.000 

472.068 

358.931 

3 

374.000 

364.345 

373.749 

364.338 

4 

477.750 

506.250 

477.665 

506.364 

5 

327.833 

507.571 

327.716 

507.587 

6 

388.000 

606.738 

387.969 

606.672 

7 

479.958 

682.204 

479.854 

682.244 

8 

269.861 

680.139 

269.660 

680.114 

9 

403.000 

842.329 

403.149 

842.443 

10 

296.000 

813.905 

296.116 

813.729 

11 

185.056 

765.151 

185.312 

765.781 

12 

63.179 

860.000 

63.179 

859.858 

13 

383.867 

297.233 

383.903 

297.258 

14 

338.917 

399.250 

338.986 

399.362 

15 

321.611 

448.611 

321.843 

448.401 

16 

266.833 

590.445 

266.906 

590.426 

17 

158.933 

720.467 

158.887 

720.118 

18 

170.681 

862.755 

170.500 

862.758 
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Table 11: Camera parameters for the Boeing dataset (example 2). Note that 
the RMS errors are given in pixels, while the orientation angles are given in 
degrees. 


Parameter 

Station 1 

Station 2 

Station 3 


2056.68 

2881.36 

2773.12 

Y c 

209.00 

201.18 

197.04 

Z c 

1961.96 

1763.88 

1814.32 

U) 

-8.821 

-9.476 

-9.434 

\ <t> 

8.948 

-18.402 

-14.592 

K 

88.806 

93.089 

-24.335 

\ ky 

0.163 

0.163 

0.163 

k 2 

-0.435 

-0.435 

-0.435 

^3 

0.901 

0.901 

0.901 

Pi 

9.063 E-3 

9.063 E-3 

9.063 E-3 

P2 

0.003 

0.003 

0.003 

a\ 

-5.496 E-3 

-5.496 E-3 

-5.496 E-3 

a 2 

4.362 E-4 

4.362 E-4 

4.362 E-4 

Xp 

-3.814 E-3 

-3.814 E-3 

-3.814 E-3 

Vv 

0.031 

0.031 

0.031 

f 

0.955 

0.955 

0.955 

x RMS Error 

0.124 

0.160 

0.143 

y RMS Error 

0.191 

0.187 

0.192 




